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The time evolution of galaxy cluster abundance is used to constrain cosmological parameters 
in dark matter models containing a fraction of hot particles (massive neutrino). We test 
the modified MDM models with cosmic gravitational waves which are in agreement with 
observational data at z = 0, and show that they do not pass the cluster evolution test and 
therefore should be ruled out. The models with a non-zero cosmological constant are in better 
agreement with the evolution test. We estimate Qa and find that it is strongly affected by a 
small fraction of hot dark matter: 0.4 < J1a < 0.8 for Qh/Qm < 0.2. 



1 Introduction 

Any realistic cosmological model should be consistent with observation data on large scale 
structure (LSS) in the Universe in the range of scales from galaxy IH (~ l/i -1 Mpc) up to large 
scale CMB anisotropy (CMBA)i (~ lOOO/i" 1 Mpc). 

Nowadays, the most popular cosmological models are cold dark matter with non-zero cosmo- 
logical constant (ACDM)@ and mixed dark matter without □ and withH cosmological constant, 
where MDM is in the form of non-baryonic hot and cold collisionless particles. 

All cosmological models were re-addressed after CMBA experimental detection! to reveal 
their positive and negative features. Both ACDM and MDM models have met several difficulties. 

As far as ACDM models are concerned i'Bi they demand a high value of the cosmological 
constant (Oa > 0.7). In this case ACDM is able to fit a set of LSS observational constraints 
whereas at small scales it overproduces the number of collapsed objects by a factor 2 in compar- 
ison with the corresponding number of gravitationally bounded objects in galaxies catalogesB. 

Regarding standard MDM models the difficulties are related vvjijj late galaxy(quasar) for- 
mation!!^ and too high number density of galaxy clusters at z = IUtEj: standard MDM is ruled 
out at 2a CL. 

One of the possibilities to overcome the sMDM difficulties is a consideration of MDM with 
some amount of cosmic gravitational wave s ( CGW) O'EJilil. The importance of fundamental 
CGW has been discussed earlier in papers Although ACDM models with < 0.3 

normalized by COBE 4-year datalia are consistent with the cluster number density test, in order 
to archive an agreement with CMBA and cluster abundance data in MDM models (for both blue 
(n > 1) and red (n < 1) scalar perturbation spectra) it is necessary to take into account the 
CGW contribution in the derived value of AT/T at 10° angular scale. The latter is estimated 
by parameter T/S, the ratio of the tensor to the scalar mode contributions. 



a n is the slope-index of post-inflationary density perturbation power spectrum. 



As an alternative to ACDM and MDM. the MDM models with a non-zero cosmological 
constant (AMDM) have been proposed intloCil. The advantage of these models is in retaining 
the inflationary paradigm and flat Harrison-Zel'dovich spectrum (n = 1) with a smaller value 
(comparing with ACDM) of cosmological A-term. Other notable features of AMDM are the 
possibility of negligible CGW contribution (T/S =0) and a small fraction of hot particles: even 
10% of massive neutrinos (Uh/Um ~ 0.1) could change the value of cosmological constant, 
being however in good agreement with other independent tests (CMB datac3, QSO lensingca, 
Hipparcos datao, SNIacl). 

In this paper we consider MDM (with CGW) and AMDM (with T/S=0), applying to these 
models the cluster evoluMnntest. The latter test constraints effectively parameter Qm (and 
hence Q,\ in flat models) Q'E§E3. Our aim is to demonstrate how the presence of hot component 
could influence the £l\ limits. 

We describe our models in Section 2 and the evolution test in Section 3, the results are 
summarized in Section 4. 



2 Cosmological Models 

We assume that DM is given by mixture of CDM and HDM components in the flat background 
space. The free model parameters are: 

• Oji, the total matter density in the Universe (fij^ = 1 — £l\ = f2f, + £Ih + f^, the latter 
are density parameters of baryons, hot, and cold particles, respectively, 

• Qh/Qm, the fraction of hot DM, 

• Nh, the number of massive neutrino species, 

• h , the Hubble constant in units 100 km s Mpc . 

The massive and massless neutrinos are described by the corresponding distribution func- 
tions which are evaluated from the Boltzmann-Vlasov collisionless equations. Cold particles 
(neutralino or hypotethical axions) are described as a pressureless fluid (pc = 0). Baryons and 
photons are treated as an ideal hydrodynamic fluid satisfying the Euler equations of motion (we 
choose the fixed value for baryon density parameter, f2& = 0.015//i 2 ). All components interact 
with each other only gravitationally. 

Assuming the power low post-inflationary density perturbations spectrum (Po(k)ook n ), the 
final total density power spectrum can be expressed as: 



P(k,z) = Ak n T 2 ' 



g(Q M (z)) 
(l + z)g(Q 



M, 



(1) 



where A is the normalization constant, T(k^z) is the total transfer function, and <?(Om(<z)) is 
the suppression coefficient c3. According toc3 g(Q,M{z)) can be approximated as 

9(!iM( ,)) = 2.5o M( ,) (4 + - sijga + 4w) <2) 

where the current matter abundance Qm{z) can be written as follows: 

n M { z ) = n M - {1 + z) * , n M = n M (o). (3) 

1 — Mm + (i + z)°Um 
We use the transfer function approximations for sCDM0 and AMDM0 models. 



3 Cluster Evolution Test 



The mass function for the gravitationally bounded halos of mass greater than M formed in the 
flat Universe by redshift z is given byH 

oo 

N(> M,z) = [ n(M',z)dM' (4) 



M 



where n(M, z)dM is the comoving number density of collapsed objects with masses lying in the 
interval (M, M + dM): 

2p5 c 1 .da(R,z) -- 



M = fpi? 3 , p is the mean matter density, and 5 C is the critical density contrast for a linear 
overdensity able to collaps. The rms amplitude of density fluctuation in the spheres of radius R 
at redshift z is related to the power spectrum as 

oo 

a\R,z) = 1^2 J P(k,z)\W(k,R)\ 2 k 2 dk, (6) 
o 

where W(kR) is a Fourier component of the top-hat window function: W(x) = (sin x — xcosx) . 

For the matter dominated Universe (Qm = 1) 5 C = 1.686 (e.g. 01). For the flat models with 
A-term S c depends weakly on the current matter abundance. The theoretical mass functions 
obtained with the help of the Press-Schechter formalism (eqs. (4,5,6)) are in a good agreement 
with other methods including numerical simulations E3a. Due to the exponential dependence 
(see eq.(5)) the cluster number density is very sensitive to the value of ag ■ 

Rich galaxy clusters are strong X-rays sources characterised by the gas temperature (see 
the review papers at this Conference). With help of numerical simulation it was confirmed that 
TpOoM 2 / 3 , the proportionality coefficient depends slightlv on cosmological models. Here, we 
make use the T — M relation for isothermal gas given in 



2 „ 1 . 1 

3 



where /?(~ 1) is the ratio of the galaxy kinetic energy to the gas thermal energy, X(~ 0.76) is 
the hydrogen mass fraction. The value A cr is the ratio of a mean halo density (within the virial 
radius of collapsed object) to the critical density of the Universe at the corresponding redshift. 
For M < 1, A cr can be derived analytically and approximated as A cr — 17809^- . 

As well as the cluster mass function A^(> M) = A^(> M, 0), the X-cluster temperature 
function A^(> T) = N(> TAT) is also sensitive to 07? = <r(R,0). The comparison of the ob- 
served temperature functions with cluster mass function Ej shows an agreement between two 
approaches. Confronting with observational data both tests provide a powerful constraint on 



Ml 



cr-value in different DM models 

The cluster mass(temperature) functions and evolution tests in ACDM has been discussed in 
details (e.g.§@@). It has been shown that flat cosmological models are preferable in comparison 
with open ones. To archive a better consistency with cluster evolution observations the value of 
Qm should be lower than that obtained from the cluster number density test at z = 0. 

The dependence of the cluster mass(temperature) function on z is model dependent due to 
the differences in growth rate of density perturbations. This is reflected in a time evolution 
related to the suppression coefficient g(p,M(z))' 

a{R,z)=a R — — — , (8) 
g{U M ) 1 + z 



where g(£l(z)) is given by eq.@. 

We consider below the cosmological models containing massive neutrino. In this case z- 
dependence also appears in the transfer function. 



4 Results and Discussion 



4.1 MDM Models with Non-Zero Tensor Mode 

To test MDM models with zero cosmological constant we follow the normalization procedure 
given in MM All models are^aonalized by ag by the best fit of observational present day cluster 
mass l 36 H 26 l and temperature E3e3 functions. The agreement with AT/T data leads to a non-zero 
derived parameter T/S. The required value of T/S is mainly determined by spectral index n, 
abundance of hot matter component , and Hubble constant h. 

Fig. la presents the present day cluster temperature functions iV(> T) for different values of 
Qh in MDM models with CGW normalized by ag ~ 0.52. The normalization does not depend 
on £Ih, therefore all curves cross each other at some fixed point (corresponding to the mass M 
in the sphere of radius 8/t _1 Mpc). 

A significant contribution of both CGW and massive neutrino is needed to fit the CMBA 
data. For h = 0.6 and £Ih = 0.2 models with flat spectra n = 1, T/S=1.7; for blue spectra 
n = 1.1, T/S=3.6; for red spectra n = 0.9, T/S=O.6 0. An increase of Q,h decreases the 
parameter T/S, however the problem arises with small-scale clustering and Ly a cloud formation 
test El. At low tin < 0.2, neither possible changes of h nor models with three species of massive 
neutrino can suppress high contribution of CGW to CMBA. Even neglecting the observational 
problem at high > 0.2 cannot help to achieve an agreement with cluster evolution test in 
MDM matter dominated models. 

Figs, lb, 2a, 2b present the cluster evolution in MDM models. None of the considered models 
can fit the data at high z > 0.3. We varied (Fig. lb.), h (Fig. 2a.), and post-inflationary 
spectral index n (Fig. 2b). All the models predict the number of galaxy clusters at high redshifts 
at least two orders of magnitude smaller than the observed one. This fact indicates that the 
evolution should be slower than that found in MDM models without a cosmological constant. 



4.2 AMD M Models without CGW 

Let us consider flat AMDM models without CGW normalized by the COBE 4-year data. 

The best fit a% obtained from nearby cluster observational abundance depends on Qm as: 
o-g(cl) = 0.52n^- 52+0A3nM It The coincidence of <r 8 (cl) with cr 8 (cmb) derived in COBE 
normalized ACDM model would limit the parameters Qm (^a) by the corresponding values 
~ 0.3(0.7) (e.g. Ell). Cluster evolution test bounds these parameters ever stronger 0, ~ 0.2(0.8). 

As an example, we show the cluster evolution functions in ACDM model with = 0.7 
for the cluster masses M > 8 • 10 14 M Q in the comoving radius 1.5h x Mpc, as functions of h 
(Fig. 3a) and n (Fig. 3b). The cluster evolution is not practically affected by changes of the 
spectral index. Regarding parameter h, the models with h > 0.65 are preferable. 

While the requested value of cosmological constant in ACDM models is quite high (Q\ > 0.7), 
a small amount of hot particles [VLh/^m < 0.2) in AMDM models allows a consistency with 
observations for smaller Q\ !0 

We performed the calculations of function A^(> T) for different parameters £l\ and Qh/Qm 
(Figs. 4a,b). Increasing Qh/^m reduces the needed value of cosmological constant. Taking into 
account the Ly a forest test, models with Qh/^m —0.1 are preferable Parameter Qm remains 
in the range 0.4 < Q,m < 0.6 in agreement with other results (cf. Bt^'EJB) . 

Figs. 5-6 present the cluster evolution in AMDM models as functions of the parameters f?A, 
Qh/&m, and h. Increasing £Ih/&m raises the needed value of J7a (in agreement with the cluster 



number density test). Changing h does not significally influence the bounds. 

Our investigation shows that cluster evolution test clearly indicate the existence of non- 
zero cosmological constant for a set of spatially flat AMDM models with negligible abundance 
of CGW. The presence of a small amount of hot particles like massive neutrino reduces the 
required value of the cosmological constant. For example, if massive neutrinos consistute only 
~ 10% of the total DM density, the models with Q A ~ 0.5 ± 0.1 and h ~ 0.65 ± 0.05 satisfy the 
observational data best. 
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Figure 1: (a) The present day cluster abundance N(> T) (data from Henry and Arnaud 1991), and (b) the 
cluster evolution N(> M = 8 ■ lO 14 M ,z) (data from Bahcall and Fan 1998), in MDM models with CGW, 
normalized as as = 0.52, with n = 1.0, h = 0.6, Q b = 0.015//i 2 , for Q H = 0,0.1,0.2,0.3 (solid, dot, short-dash, 

long dash lines, resp.). The needed T/S is shown. 




Figure 2: The cluster evolution N(> M = 8- 1O 14 M , z) in MDM models normalized as as = 0.52 with Qh = 0.2, 
fib = 0.015/ft 2 , for (a) n = 1, h = 0.5, 0.6, 0.7 (solid, dot, short-dash lines, resp.), and (b) h = 0.6, n = 0.9, 1.0, 1.1 
(solid, dot, short-dash lines, resp.). The needed T/S is shown. The points correspond to Bahcall and Fan 1998. 
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Figure 3: The cluster evolution N(> M = 8-lO 14 M , z) in ACDM models normalized as cr 8 = 0.52fi^ - 52nM+013) 
with Qb = 0.015/ft 2 , Qa = 0.7, for (a) n — 1 and ft = 0.65,0.70,0.75 (solid, dot, short-dash, long-dash lines, 
resp.), and (b) h = 0.7, and n — 0.8,0.9, 1.0, 1.1 (solid, dot, short-dash, long-dash lines, resp.). The data points 

correspond to Bahcall and Fan 1998. 
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Figure 4: The present day cluster abundance iV(> T) in AMDM models normalized by COBE 4-year data with 
h = 0.6, Q b = 0.015//I 2 , n= 1, for tt A = 0,0.31,0.45,0.65,0.74 (solid, dot, short-dash, long-dash, dot-dash lines, 
resp.), and (a) Hh/Qm = 0.1, (b) Qh/Qm = 0.2. The points correspond to Henry and Arnaud 1991. 




Figure 5: The cluster evolution N(> M = 8-lO 14 M , z) in AMDM models normalized as a s = 0.52^7 u ' 52SiM+u '" ) 
with Qb — 0.015/fo 2 , h = 0.6, n = 1, for Qa = 0,0.31,0.45,0.55,0.65 (solid, dot, short-dash, dot-dash lines, long- 
dash, resp.), and (a) Qh /CIm = 0.1, (b) Qh/£Im = 0.2. The data points correspond to Bahcall and Fan 1998. 




Figure 6: The cluster evolution JV(> M = 8 ■ 10 14 M Q ,z) in AMDM models normalized as erg = 0.52n^ - 52+012n 
with Ov = 0.55, Q b = 0.015//i 2 , n = 1, for (a) h = 0.6, Vt H /0. M = 0,0.1,0.2 (solid, dot, short-dash, resp), and 
(b) Qh/Qai = 0.1, h = 0.6,0.65,0.7 (dot, short dash and long dash lines, resp.). The data points correspond to 

Bahcall and Fan 1998. 
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